Searches for New Physics in Diphoton Events in pp collisions at yfs =1.8 TeV 
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We have searched for anomalous production of missing Et (Et)j jets, leptons (e, /i, r), 6-quarks, or 
additional photons in events containing two isolated, central < 1.0) photons with Et > 12 GeV. 
The results are consistent with standard model expectations, with the possible exception of one 
event that has in addition to the two photons a central electron, a high-Ex electromagnetic cluster, 
and large E T . We set limits using two specific SUSY scenarios for production of diphoton events 
with E T . 



PACS numbers 13.85Rm, 13.85Qk, 14.80.-j,14.80.Ly 
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In many models involving physics beyond the standard model (SM), cascade decays of heavy new particles generate 
77 signatures involving missing transverse energy (^1 T ), jets, leptons, gauge bosons (W, Z , 7), and possibly b- 
quarks jjj. For example, in supersymmetric models with a light gravitino, pair-production of selectrons which decay 
via e — ► eNi — > e-fG produces the 77 final state along with E T and electrons. In the data taken during 1993-1995, 
an 'ee77^ T ' candidate event was recorded with the CDF Detector ||. We have performed a systematic search 
for other anomalous 77 events by examining events with two isolated, central (|ry| < 1.0) photons with Et > 12 GeV 
which contain E T , jets, leptons (e,n, r), 6-quarks, or additional photons 0. This search is based on 85±7 pb _1 of 
data from pp collisions at y/s = 1.8 TeV collected with the CDF detector. In this Letter we describe the results of 
the search, including the ee77E T candidate event, and set limits on two SUSY models that have arisen to explain it. 

We briefly describe here the relevant aspects of the CDF detector. The magnetic spectrometer consists of tracking 
devices inside a 3-m diameter, 5-m long superconducting solenoidal magnet which operates at 1.4 T. A four-layer 
silicon microstrip vertex detector (SVX) 0] is used to identify b hadron decays. A set of vertex time projection 
chambers (VTX) surrounding the SVX is used to find the z position of the pp interaction (z vertex ). The 3. 5-m long 
central tracking chamber (CTC) is used to measure the momenta of charged particles. The calorimeter, constructed 
of projective electromagnetic and hadronic towers, is divided into a central barrel which surrounds the solenoid 
coil (|^| < 1.1), 'end-plugs' (1.1 < \r/\ < 2.4), and forward/backward modules (2.4 < \rj\ < 4.2). Wire chambers with 
cathode strip readout give 2-dimensional profiles of electromagnetic showers in the central and plug regions (CES and 
PES systems, respectively). A system of proportional wire chambers (CPR) in front of the central electromagnetic 
calorimeters uses the 1-radiation- length-thick magnet coil as a 'preradiator', allowing photon/7r° discrimination on a 
statistical basis by measuring the conversion probability Muons are identified with the central muon chambers, 
situated outside the calorimeters in the region \r)\ < 1.1. 

The data sample selection starts with events with two photon candidates identified by the three- level trigger [Q. 
At Level 1, events are required to have two electromagnetic calorimeter trigger-towers || with measured Et of more 
than 4 GeV. At Level 2, we require the logical 'OR' of two triggers, one optimized for good background rejection at 
low Et and the other for high efficiency at high Et- The low-threshold diphoton trigger requires two electromagnetic 
clusters || with Et > 10 GeV and an isolation requirement of less than 4 GeV in a 3-by-3 array of trigger-towers 
around the cluster; the high-threshold (16 GeV) trigger has no isolation requirement. Corresponding Level 3 triggers 
require cluster energies calculated with the offline photon algorithm Q to be above the 10 GeV and 16 GeV thresholds. 
The low-threshold trigger also requires the clusters be in a restricted fiducial region of the calorimeter |l(| . 

We use the following selection criteria offline: a) two isolated [[ll| central electromagnetic clusters with Et> 12 GeV 
(where the 10 GeV trigger becomes > 98% efficient); b) no tracks, or only one track with Pt< 1 GeV, pointing at 
either cluster (to remove electrons or jets); c) pulse height and shape in the CES consistent with a shower due to 
a photon (to remove ir° backgrounds); d) no other photon candidate within the same 15° segment of the CES (to 
remove n° backgrounds) jf2j; e) |z ver t cx | < 60 cm (to maintain the projective geometry of the calorimeter); and f) no 
energy out-of-time with the collision (to suppress cosmic rays) |13j . For events in which both photon candidates 
have Et> 22 GeV (where the 16 GeV trigger becomes > 98% efficient) the fiducial and isolation requirements are 
relaxed |lO|]llfl . The final data set consists of 2239 events. 

The efficiency for identifying an isolated photon is measured using electrons in a control sample of 2663 Z° — > e + e~ 
decays to be 68 ± 3% for the 12 GeV selection criteria, and 84 ± 4% for the 22 GeV selection criteria, in each 
case approximately flat in Et- Photon backgrounds are measured using the shower shape in the CES system for 
Et< 35 GeV, where the difference between a single 7 or 7r° — * 77 can be resolved, and the conversion probability in 
the CPR for Et> 35 GeV ||. Since the purity of the 77 sample is of less importance to searches for rare signatures 
than the efficiency, we have chosen selection criteria to keep the efficiency high; however these admit a substantial 
number of background events. The fraction of events in the sample which contain two prompt photons is measured 
to be 15 ± 4%. 

We search the diphoton events for the presence of E T , jets, electrons, muons, taus, 6-quarks, and additional 
photons. To minimize the effect of fluctuations due to mismeasurement of jet energies, we recalculate the E T by 
making jet energy corrections which take into account cracks between detector components and nonlinear calorimeter 
response |l4||ll|. In the 77 + E T search we remove events which have a jet with uncorrected Et> 10 GeV pointing 
within 10 degrees in azimuth of the E T . The resulting resolution on either the x or y component of E T , determined 
from a study of Z° bosons, is well-parameterized by cr(E T ) = (2.66 ± 0.34 GeV) + (0.043 ± 0.007) x EE T , where SE T 
does not include the Et of either electron. The criteria for identifying jets (uncorrected Et> 10 GeV and < 2.0), 
electrons, muons and 6-jets are identical to those used in the top-quark discovery [T^ ]. The tau selection is the same 
as used in the study of it decays into er and \xt final states Jl(|. Any third photon is required to have Et> 25 GeV 
and to pass the high-threshold selection criteria. 

Table H summarizes the observed and expected numbers of events. The distributions in E T and the number of 
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jets, Njet, are shown in Figure [|. The shapes of the E T distributions are in good agreement with the resolution 
derived from the Z° control sample, shown as the hatched region in Figures pi and [l]c. The distributions in Nj et are 
well-modeled by an exponential extrapolation, shown in Figures |b and 0d. For a photon threshold of 12 GeV we 
observe 1 event with E T > 35 GeV, with a SM expectation of 0.5 ± 0.1, and 2 events with 4 or more jets, versus an 
expectation of 1.6 ± 0.4. For a photon threshold of 25 GeV, we observe 2 events with E T > 25 GeV, with 0.5 ± 0.1 
expected, and events with 3 or more jets, with 1.7 ± 1.5 expected. 

We find 2 events with 6-tags, consistent with background expectations, and no events with a third photon. There 
are 4 events with a central lepton: one event is consistent with a double-radiative Z° decay (m MM77 = 92 ± 1 GeV/c 2 ), 
one is consistent with a radiative Z° decay with a lost track (m ee7 = 91 ± 2 GeV/c 2 ), and one has a r candidate, for 
which we expect a fake background of 0.2 ± 0.1 events. From Table || and Figure [l], we find agreement between our 
observations and SM model predictions with one possible exception. 

The event that has the largest E T (E T — 55 ± 7 GeV) among all diphoton candidates, has in addition to the two 
high-Ex photons a central electron and an electromagnetic cluster in the plug calorimeter which passes the electron 
selection criteria used for Z° identification JT^] . The 4- vectors are presented in Table 0. Because the momenta of the 
four clusters are measured by the electromagnetic calorimeters, the resolution on each is a few per cent (see Table ||). 
The total Pt of the 4-cluster system is 48 ± 2 GeV/c, opposite to the E T and in good agreement with the measured 
magnitude, implying the imbalance is intrinsic to the 4-cluster system. The invariant mass of the electron and the 
electromagnetic cluster in the plug calorimeter is 163 ± 3 GeV/c 2 , far from the Z° mass. The invariant mass of the 
4-body system is 232 ± 4 GeV/c 2 ; a lower limit on the invariant mass of the total system is found by including the 
E T (taking Pz = 0), to be 307 ± 9 GeV/c 2 . 

Although the electromagnetic cluster in the plug calorimeter passes all of the standard electron selection criteria jl7) , 
there is no track in the SVX pointing directly at the cluster, as would be expected if the cluster were due to an 
electron ^||l^]. There is, however, a track 26 mrad away in <fi. Using a sample of 1009 electrons in the end-plug 
calorimeter from Z {) — > e + e~ events, we estimate the resolution on <f> to be 1.5 mrad; no events have a mismatch 
of greater than 20 mrad. The probability of an electron to have a <fi mismatch this large is thus less than 0.3% at 
95% C.L. The interpretation of the cluster as coming from an isolated photon, the 1-prong hadronic decay of a t, or 
a jet, while possible, are also all unlikely in that this would be an unusual example of any of them [Q. We simply do 
not have enough information to establish the origin of the cluster. 

We have estimated the SM rates for producing a signature of two photons, two electromagnetic clusters (one central) 
passing the electron requirements and E T , all with Et > 25 GeV, and m ee > 110 GeV/c 2 (above the Z° boson) [Q. 
Using both data and Monte Carlo methods, we have considered production of SM WW'yy and ti, as well as sources 
which include additional cosmic ray interactions, jets which fake electrons and/or photons, and overlapping events. 
The total rate is 1 x 10~ 6 events, with the dominant sources being WW-f-f (8 x 10~ 7 events) and WWjj (8 x 10~ 8 
events). Removing sources where the plug cluster is due to a real electron, the rate is reduced to 6 x 10~ 8 events, with 
the dominant source being Wfyj (5 x 10~ 8 events). Multiple events in the same beam crossing |l{| or the overlap 
of a cosmic ray interaction with a pp event contribute a total of 8 x 10~ 9 events. We emphasize that while these SM 
estimates are small and have led to valuable speculation M, it is indefensible to claim evidence of new physics based 
on one peculiar event selected from 3 x 10 12 events. 

One possible source of an ee77E T signature is anomalous VFW77 production. This hypothesis can be checked 
by searching for events where each W decayed hadronically rather than leptonically. A Monte Carlo study using a 
standard model WWjj calculation |2(J shows that anomalous WW77 production would produce detected events 
with two photons with Et>25 GeV and three or more jets 30 times more often than events with two photons, two 
leptons and E T . No events with three or more jets are seen in the Nj et distribution (Figure [j]d). 

We proceed to set limits on two SUSY models. There has been recent interest in supersymmetric models with either 
the lightest neutralino (Ni) decaying into a photon and gravitino (G) [16, lc, le, lg], N\— > 7G, or a supergravity 
scenario in which the second- lightest neutralino decays via a loop into the lightest neutralino and a photon Q]c, 
./V2— ► 7V1. Both of these models would produce events with two photons and E T . 

We use the SPYTHIA Monte Carlo |^l[ with a full detector simulation to investigate the N2 — ► jNi model of 
Ambrosanio et al. with Mjy = 36.6 GeV and M^y 2 = 64.6 GeV f22fl . Direct production and cascade decays are 
predicted to produce 2.4 events that pass the selection criteria of E^ > 12 GeV and E T > 35 GeV. In the data only 
the one event passes this selection; we consequently cannot exclude this model. To provide a normalization point 
{e.g. for model-builders to estimate the detector efficiency), we have simulated direct N2N2 production for this same 
model and find an acceptance of 5.4%. Treating the one event as signal, and performing no background subtraction, 
we derive a 95% C.L. cross section upper limit of 1.1 pb. 

Production of 77 events in the light gravitino scenario of Babu et al. We is dominated by C1N2 and C\C\ production 
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and decay. Figure g shows the cross-section limits, using the same methods, versus the mass of the C\. The lines 
show the experimental limit and the theoretically predicted cross section for the lowest value of that is excluded 
(Mcj < 120 GeV at 95% C.L., for tan/3 = 5, fx < 0). Note that because the C\ and N\ masses are related we also 
exclude Mjvi < 65 GeV at 95% C.L. {n > 0, tan/3 = 5). These limits are similar to those of the D0 collaboration p3| . 

In conclusion, we have searched a sample of 85 pb _1 for events with two central photons and anomalous production 
of missing transverse energy, jets, charged leptons (e, fi, and r), 6-quarks and photons. We find good agreement with 
standard model expectations, with the possible exception of one event which has unusually large E T and in addition 
to the two photons has a high-ET central electron and a high-Ex electromagnetic cluster. 
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Education of Japan, the Natural Sciences and Engineering Research Council of Canada, the National Science Council 
of the Republic of China; and the A. P. Sloan Foundation. 
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Signature (Object) 


Obs. 


Expected 


Ref. 


E T >35 GeV, |A0p T _ jot | > 10° 


1 


0.5 ±0.1 






Njct > 4, > 10 GeV, |j7 jot | < 2.0 


2 


1.6 ± 0.4 


I 


jj 


6-tag, > 25 GeV 


2 


1.3 ±0.7 




jj 


Central 7, EiJ 3 > 25 GeV 





0.1 ±0.1 






Central e or /1, E; or " > 25 GeV 


3 


0.3 ±0.1 


I 




Central r, E T > 25 GeV 


1 


0.2 ± 0.1 




m 


TABLE I. Number of observed and expected 77 events with additional objects in 85 pb . The selection criteria, efficiencies, 


and background estimation methods used 


in identifying the jets, leptons and fo-tags are discussed in the references. 






Px 


Py 


P z E 


Et 


(GeV/c) 


(GeV/c) 


(GeV/c) (GeV) 


(GeV) 


71 32.1(9) 


-16.8(5) 


-35(1) 50(1) 


36(1) 


72 -12.9(4) 


-29.6(9) 


-22.5(7) 39(1) 


32.3(9) 


e" -34(1) 


11.5(3) 


21.7(6) 42(1) 


36(1) 


Plug EM Cluster 60(2) 


19.0(5) 


-172(5) 183(5) 


63(2) 


E T -54(7) 


13(7) 




55(7) 



TABLE II. The 4-vectors of the objects in the eeyyf) T candidate event. The parentheses represent the uncertainty in the 
last digit. There are no jets with Et > 10 GeV. 
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FIG. 1. a) The E T spectrum for events with two central photons with Et> 12 GeV. We have removed events which have 
any jet with Et> 10 GeV pointing within 10 degrees in azimuth of the E T . The cross-hatched regions represent the background 
estimates derived from the E T resolution in the Z° control sample, b) The spectrum in number of jets with Et> 10 GeV and 
\rj\ < 2.0 (Njot) for events with two central photons with Et> 12 GeV. c) and d) the same plots with photon Et> 25 GeV. 
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FIG. 2. The cross-section upper limits versus the mass of the C\ for the light gravitino scenario of Babu et al. jjje. The shaded 
region shows the range of cross section limits as the parameters are varied within the ranges 1 < tan/3 < 25, M2 < 200 GeV, 
and n > or fj, < 0. The lines show the experimental limit and the theoretically predicted cross section for the lowest value of 
M Cl that is excluded (M Cl < 120 GeV at 95% C.L., for tan/3 = 5, n < 0). 
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